#!/usr/bin/env python
import catsHTM
import os
import scipy as S
import sys
#-
#- chao wu<cwu@bao.ac.cn>
#- to search ps1.
#-

ra=float(sys.argv[1])  #- deg
dec=float(sys.argv[2]) #- deg
radius=float(sys.argv[3]) #- arcsec


ra1=  S.deg2rad(ra)
dec1= S.deg2rad(dec)
dd= catsHTM.cone_search('PS1',ra1,dec1,radius,catalogs_dir=os.getenv('CHTMPATH'),OnlyCone=False)
dd0=dd[0].copy()
dd1=dd[1].copy()
dd2=dd[2].copy()

dd0[:,0]=S.rad2deg(dd[0][:,0])
dd0[:,1]=S.rad2deg(dd[0][:,1])

dd2[0]='deg'
dd2[1]='deg'

tt=" ".join(list(dd1))
print ("#- ra=%s dec=%s radius(arcsec)=%s " %(ra,dec,radius))
print ("#-%s" %tt)
for d0 in dd0:
    d0s = " ".join(list(d0.astype('S10')))
    print (d0s)

#- output format:    
##- array(['RA', 'Dec', 'ErrRA', 'ErrDec', 'MeanEpoch', 'posMeanChi2',
#        'gPSFMag', 'gPSFMagErr', 'gpsfLikelihood', 'gMeanPSFMagStd',
#        'gMeanPSFMagNpt', 'gMeanPSFMagMin', 'gMeanPSFMagMax', 'rPSFMag',
#        'rPSFMagErr', 'rpsfLikelihood', 'rMeanPSFMagStd', 'rMeanPSFMagNpt',
#        'rMeanPSFMagMin', 'rMeanPSFMagMax', 'iPSFMag', 'iPSFMagErr',
#        'ipsfLikelihood', 'iMeanPSFMagStd', 'iMeanPSFMagNpt',
#        'iMeanPSFMagMin', 'iMeanPSFMagMax', 'zPSFMag', 'zPSFMagErr',
#        'zpsfLikelihood', 'zMeanPSFMagStd', 'zMeanPSFMagNpt',
#        'zMeanPSFMagMin', 'zMeanPSFMagMax', 'yPSFMag', 'yPSFMagErr',
#        'ypsfLikelihood', 'yMeanPSFMagStd', 'yMeanPSFMagNpt',
#        'yMeanPSFMagMin', 'yMeanPSFMagMax']
        
